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ABSTRACT 

We analyze predictions from two independently developed galaxy formation models to 
study the mechanisms, environments, and characteristic times of bulge formation in a 
ACDM cosmogony. For each model, we test different prescriptions for bulge formation 
in order to quantify the relative importance of different channels. Our results show 
that the strong correlation between galaxy and halo mass for central galaxies, and the 
richer merger history of more massive systems naturally give rise to a strong correla- 
tion between galaxy mass and morphology, and between halo mass and morphological 
type of central galaxies. Large fractions of the bulge mass are acquired through major 
and minor mergers, but disc instability plays an important role, particularly for inter- 
mediate mass galaxies. We find that the modelling of disc instability events, as well as 
of the galaxy merger times, can affect significantly the timing of bulge formation, and 
the relative importance of different channels. Bulge dominated galaxies acquire their 
morphology through major mergers, but this can be modified by cooling of gas from 
the surrounding hot halo. We find that disc regrowth is a non negligible component of 
the evolution of bulge dominated galaxies, particularly for low to intermediate masses, 
and at high redshifts. 

Key words: galaxies: formation - galaxies: evolution - galaxies: bulges - galaxies: 
interactions - galaxies: structure. 



1 INTRODUCTION 

Perhaps the clearest and most convenient definition of a 
bulge is that of a centrally concentrated stellar distribu- 
tion, with a smooth and spherical appearance. Indeed, such 
a definition unde rlies the classification scheme introduced 
bv lHubblei { 192rj ). In the local Universe, about 60 per cent 
of the total stellar ma ss of massive g alaxies is contained in 
ellipticals and bulges jGadotti 2009). It is clear then, that 
understanding how bulges form and evolve is integral to the 
question of understanding galaxy formation and evolution. 

Until the early 1980s, bulges were thought to belong to 
the same family as elliptical galaxies, and to have formed 
through the same physical process(es). Several lines of evi- 
dence, however, indicated later that the class 'bulges' rep- 
resents a heterogeneous family including systems with very 
different properties, and likely very different formation and 
evolutionary histories. Indeed, significant differences were 
found between the kinematics of ellipticals and bulges (e.g. 
iDressler fc Sandagej 1 19831 : iDavies fc Illingworthl Il983l . and 
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references therein). It was also noted that many bulges ex- 
hibit a 'boxy' or 'peanut' shaped structure at small radii. 
This shape, that is unlikely to be due to the gravita- 
tional influence of the disc, was fo und to be associated with 
differ ential cylindrical rotation (|Kormendv fc Illingworthl 
1982). In the past decades, substantial evidence has ac- 
cumulated that many bulges have 'disc-l i ke' exponen- 
tial p r ofiles (e.g. |Andredakis fc Sandersll 19941; ICarollo et al.l 



l200ll : iBalcells et al l 120031 ; I Fisher fc Droryl l200ct) and, in 
some cases, 'disc-like' cold kinematics (iKormendvl Il993l : 
iPinknev et"al]|2003l ; Ikormendv fc Kennicuttll2004 and ref- 
erences therein). These systems are now often referred to 
as 'pseudo- bulges', as opposed to 'classical bulges' that are 
relatively featureless both dynamically and photometrically, 
and appear to have a close affinity with elliptical galaxies. 

The current view is that classical bulges are formed 
through rapid collapse or hierarchical mergers of smaller ob- 
jects, and corresponding dissipative gas processes. Early nu- 
merical simulations showed that close interactions can lead 
to a strong internal dynamical response, driving the forma- 
tion of spiral arms and, in some cases, of strong bar modes. 
The axisymmetry of these structures induces a compression 
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of the gas that can fu el nuclear st arbursts and/or nuclear 
AGN activity (see e.g. lMihosll2004l . and references therein). 
Simulations have also shown that the merger of two spiral 
galaxies of comparable mass can produce a remnant with 
structural and photo metric properties resembling those of 
elliptical galaxies (e.g.lToomre fc Toomrdll972l : lMihosll2004l ; 
ISpringel. Pi Matteo fc Hernquisdl2005l ). On the other hand, 
pseudo-bulges are thought to originate from the evol ution 
of disc instabilities such as bars. Early simulations bv lHohll 
showed that bar formation is accompanied by a rear- 
rangement of disc material, which results in the formation of 
a high-density central core. Later and more detailed simula- 
tions have confirmed that gravitational instabilities such as 
spirals and bars are able to build 'bulge-like' structures, ei- 
ther through vertical resonances or through bending ('buck- 
ling') of the bars (e.g.lCombes et alj[l990t iRaha et al.lll99ll ; 
iDebattista et al.ll200fi ~ 

All these processes are at play in the general framework 
of hierarchical galaxy formation: galaxies are supposed to 
form through the condensation of gas at the centre of dark 
matter haloes. Conservation of angular momentum leads to 
the formation of a rotationally supported disc. If the cooling 
is 'rapid' (at high redshift and in relatively small haloes), the 
short dynamical times lead to an intense star-burst activity. 
Mergers and instabilities form 'bulges', that can eventually 
grow a new disc, provided the system is fed by an apprecia- 
ble cooling flow. In this framework then, bulge-dominated 
galaxies can be 'transitory' systems. The importance of disc 
regrowth and its correlation with the physical properties 
and/or environment of galaxies has, however, not been an- 
alyzed in detail. 

Accurate studies of the structural and physical proper- 
ties of bulges and ellipticals are now being carried out (e.g. 
lGadottill2009l . at low redshift). These studies and their ex- 
tension to higher redshift, will likely provide important con- 
straints on how the different population of bulges evolved 
as a function of cosmic time. It is therefore interesting to 
analyze in more detail predictions from recently published 
galaxy formation models, with particular regard to the ques- 
tion of what is the relative role of different physical mech- 
anisms (e.g. mergers vs disc instability) in the formation of 
galaxy bulges, and their evolution as a function of redshift, 
environment and galaxy mass. 

In this paper, we analyze predictions from two indepen- 
dently developed semi-analytic models of galaxy formation. 
While the models used in this study do not allow a fine 
classification into 'bulges', 'pseudo-bulges', or 'bars' to be 
made, they allow us to quantify the amount of mass that is 
contributed to the spheroidal components by different 'chan- 
nels' (minor and major mergers, and disc instability), and 
to study when and in which environment (s) these processes 
take place. Using two different models and, within them, dif- 
ferent prescriptions for the formation of bulges, we are able 
to analyze how the relative importance of different channels 
varies as a function of different specific physical assump- 
tions. Some of these issues have been addre ssed using similar 
classes of models in previous s tudies l|Parrv. Eke fc Frenkl 
2009; B cnson fc Devereuxll201Ch , and we will comment on 
these results below. In this study, we focus on theoretical 
predictions, and defer a detailed comparison between model 
results and observational data to a future work. 

The layout of the paper is as follows. In Section [2] we 



introduce the galaxy formation models used in this study, 
focusing on those aspects of the models that are relevant for 
bulge formation. In Section [3] we discuss the basic trends 
predicted as a function of the galaxy stellar mass and of 
the parent halo virial mass. In Sections [4] and O we analyze 
the times and environments that characterize the formation 
of galaxy bulges through different channels. In Section [6] 
we study the formation history of 'elliptical' galaxies and 
address the issue of disc regrowth. Finally, we discuss our 
results, and give our conclusions in Section [7] 



2 THE GALAXY FORMATION MODELS 

In this paper, we consider predictions from two dif- 
ferent and independently developed semi-analytic mod- 
els of galaxy formation within a ACDM cosmogony. In 
particular, we use (i) the recent impl e menta tion of the 
Munich model by |Pe Lucia fc Blaizotl (120071 ). with its 
generalization to the WMAP3 cosmology discussed in 
I Wang et afl (120081); and (ii) the MORGA NA model presented 
in iMonaco. Fontanot fc Taffonil l|2007l), and adapted t o a 
WMAP3 cosmology as described in lLo Faro et ail l|2009l ). 

Comparisons between different specific predictions from 
these two models have been di scussed in iFontanot et al.l 
(|20Qgh and lFontanot et al.l l|2010l ). A detailed analysis of the 
prescriptions adopted to model gas cooling and galaxy merg- 
ers is given in lDe Lucia et all (|2010h - We refer to the original 
papers for more details on the modelling of various physical 
processes. In the following, we highlight the main differences 
between the implementations of these ingredients, focusing 
on those physical processes that are relevant for bulge for- 
mation. We al so summarize the prescriptions proposed by 
I Hopkins et al.l (2009a), that we have implemented in the two 
models used in this study. 



2.1 The WDL08 model 

Cosmologica l fram ework: The model discussed in 
I Wang et ail l|2008l . WDL08 hereafter) takes advantage 
of N-body simulations that follow the evolution of 
N — 540 3 particles within a comoving box of size 
125/i _1 Mpc on a side. This corresponds to a parti- 
cle mass of 7.78 x 1O 8 /i _1 M0. In this study, we use 
their simulation with WMAP3 cosmological parameters: 
n m = 0.226, fi b = 0.04, fi A = 0.774, n = 0.947, and 
erg = 0.722. The Hubble constant is parametrized as 
Ho = 100 /ikms -1 Mpc -1 , and this particular simulation 
assumes h = 0.743. 

Merger trees: Simulation data were stored at 64 out- 
put times, that are approximately logarithmically spaced 
between z=20 and 1, and approximately linearly spaced in 
time thereafter. Group catalogues were constructed using a 
standard friends-of-friends (FOF) algorithm, with a linking 
length of 0.2 in units of the mean particle separation. Each 
group was then decomposed into a se t of disjoint substruc - 
tures using the algorithm SUBFIND (|Springel et al.ll200ll ). 
which iteratively determines the self-bound subunits within 
a FOF group. The most massive of these substructures is 
often referred to as the main halo, while this and all other 
substructures are all referred to as subhaloes or substruc- 
tures. Only subhaloes that retain at least 20 bound parti- 
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cles after a gravitational unbinding procedure are considered 
'genuine' subhaloes, therefore setting the subhalo detection 
limit to 2.22 x 1O 1O M0. These catalogues were then used 
to construct merger history trees of all gravi tationally self- 
bound substru ctures, as explained in deta il in lSpringel et al.l 
i|2005l . see also lDe Lucia fe Blaizotll2007l ). 

Galaxy mergers: At variance with the other model used 
in this study, the WDL08 one follows dark matter haloes af- 
ter they are accreted onto larger systems. This allows the 
dynamics of satellite galaxies residing in infalling haloes 
to be properly followed, until the parent dark matter sub- 
structure is 'destroyed' (i.e. falls below the resolution limit 
of the simulation) by tidal truncatio n and stripping (e.g. 
|Pe Lucia et alj|2004 iGao et al.ll2004) . When this happens, 
galaxies are assigned a residual surviving time using the clas- 
sical dynamical friction formulcQ. The residual merging time 
is estimated from the relative orbit of the two merging ob- 
jects, at the time of subhalo disruption. 

In the case of a 'minor' merger, the stellar mass of 
the merged galaxy is transferred to the bulge component 
of the remnant galaxy, and a fraction of the combined 
cold gas from both galaxies is turned into stars as a re- 
sult of the merger. The efficiency of the merger-driven star- 
burst is paramet r ized a dopting the formulation proposed by 
ISomerville et all l|200ll) : 

eburst = /9burst(m 2 /mi) QburBt 

where 7712/7711 is the baryonic (gas + stars) mass ratio, and 
Qburst = 0.7 and Eburst = 0.56 have been chosen to provid e 
a good fit to the numerical simulations of I Cox et al. (2008). 

All stars that form during the burst, as well as all re- 
maining cold gas, are added to the disc of the remnant 
galaxy. If the baryonic mass ratio of the merging galaxies is 
larger than 0.3, we assume that we witness a 'major' merger, 
that gives rise to a more significant starburst and destroys 
the disc of the central galaxy completely, leaving a purely 
spheroidal stellar remnant. The remnant galaxy can grow a 
new disc later on, provided it is fed by an appreciable cooling 
flow. 

Disc instability: bulges can also grow through disc insta- 
bilities, that are ass umed to take place when the followin g 
condition is verified (Efstathi ou. Lake fc Negropontell 1982ft : 

(Gm disc /r disc )V2 ~ eiim U 

In the above equation, 77idi sc , ?"disc, and Vdisc are the stel- 
lar mass, the radius, and the velocity of the disc, respec- 
tively. In this model, Vdisc = V m &x, and is computed di- 
rectly from the underlying iV-body simulation; rdi sc is the 
half-mass radius of the disc that, for an exponential disc, 
is equal to 1.68 x Rd\ Rd is the disc s cale l ength, and is 
computed following iMo. Mao fc White! (|l998l ). The model 
assumes t\i m = 0.75, that is chosen in order to reproduce 
the observed morphological mix in the local Universe. For 
each galaxy, and at each time-step, we check whether the 
instability condition is verified and, when this is the case, 
we transfer enough stellar mass from the disc to the bulge 
so as to restore stability. 

For a detailed discussion of the adopted form ulation, and for 
a com parison with different implementations, see lDe Lucia et al.l 
j201Ch . 



2.2 The MORGANA model 

Cosmological framework: The results from the MORGANA 
model presented in this study have been obtained using a 
144/i _1 Mpc box with N = 1000 3 particles, and adopting a 
cosmology with O m = 0.24, ft A = 0.76, erg = 0.8, n = 0.96, 
and h = 0.72. The dark matter data us ed by MORGANA ar e 
obtained using the code PINOCCHIO l|Monaco et al.ll2002l ). 
This algorithm, based on Lagrangian perturbation theory, 
has been shown to provide mass assembly histories of dark 
matter haloes that are in e xcellent agreem ent with results 
from numerical simulations (|Li et al.ll2007i ). 

Merger trees: For details on the con struct ion of merger 
trees, we refer to I Monaco et alj (|2007l ) and iTaffoni et al.l 
(2002). We note that PINOCCHIO does not provide infor- 
mation on dark matter substructures, so MORGANA is es- 
sentially based on the equivalent of FOF merger trees. 

Galaxy mergers: In order to model the orbital decay of 
dark matter subhaloes and galaxy mergers, MORGANA uses 
a s lightly updated vers ion of the fitting formulae provided 
by ITaffoni et al.l l|2003h . These take into account dynami- 
cal friction, mass loss by tidal stripping, tidal disruption of 
substructures, and tidal shocks. In practice, whenever two 
(FOF) haloes merge, the galaxy associated with the smaller 
halo is assigned a galaxy merger time by interpolating be- 
tween the two extreme cases of a 'live satellite' (where the 
object is subject to significant mass losses) and that of a 
'rigid' satellite (that does not suffer a significant mass loss). 
We refer to the original paper for details on the implemen- 
tation. 

As in WDL08 , MORGANA distinguishes between minor 
and major galaxy mergers, using the same baryonic mass ra- 
tio threshold (0.3). During a minor merger, the stellar mass 
and the cold gas of the accreted satellite are added to the 
bulge component of the remnant galaxy, whose disc is unaf- 
fected by the merger. During major mergers, the stellar and 
gaseous disc of the remnant galaxy are destroyed and relaxed 
into a single spheroidal component. The cold gas associated 
with the bulge can be efficiently converted into stars, and 
this occurs on very short time-scales (effectively triggering a 
'starburst') during major mergers. As in WDL08 , the rem- 
nant galaxy can grow a new disc, out of the gas cooling at 
later times. 

Disc instability: For this process, MORGANA adopts the 
same stability criterion as in the WDL08 model, but uses 
different definitions for the ma ss, radius and veloci ty of the 
disc, and assumes en m = 0.7 (|Lo Faro et al.|[2009T ). As for 
WDL08 , this is chosen in order to reproduce the observed 
morphological mix in the local Universe. In this model, mdi sc 
is the total baryonic mass of the disc, rdisc is the disc scale- 
length (also computed following IMo et alH l998). and Vdisc 
is the rotational velocit y of the disc, computed as detailed 
in iMonaco et al] l|2007i ). When the instability condition is 
verified, half of the baryonic mass of the disc is transferred 
to the bulge component. As explained above, the presence 
of a significant amount of cold gas in the bulge can trigger 
a burst of star formation. 

Additional processes: MORGANA includes additional 
physical mechanisms that influence the assembly history of 
bulges. In particular, the model allows infall of gas onto an 
existing bulge, by a fraction equal to the fraction of disc 
mass embedded in the bulge. In addition, the model takes 



© 2010 RAS, MNRAS OOO.rTHlSl 



4 G. De Lucia et al. 



into account tidal stripping of stars in satellites, and assumes 
that a fraction (f sca — 0.7 in the standard model) of satellite 
stars are unbound during major mergers and incorporated 
into a 'diffuse' light component. In the following, we will ne- 
glect the second process (i.e. we will assume f 3ca = 0). We 
have verified, however, that this does not alter significantly 
any of the results discussed in this study. 

2.3 Model differences 

The previous sections clarify that the two models used in 
this study differ in a number of details. As we will see in the 
following, some of these are reflected in significant differ- 
ences between model predictions. In this section, we briefly 
comme nt on these expe c tation s . 

In |Pe Lucia et al.l (|2010T ). we compared the merger 
model adopted in MORGANA with that employed in WDL08 
, and showed that the former provides merger times that 
are systematically shorter (by an order of magnitude) than 
those predicted by the latter. This will likely translate into 
a different relative importance of the merger channel. In or- 
der to quantify the significance of this different treatment of 
galaxy mergers, we will also show results from MORGANA 
obtained using longer merger times. 

Another significant difference between the two models 
is related to the adopted treatment of disc instability: al- 
though both models ar e based on the criterion proposed by 
Efsta thiou et al.l (| 1982T ) . they adopt different definitions for 
the mass, radius and velocity of the disk, and instability 
events have rather different consequences. In WDL08 , disc 
instability is evaluated only for the stellar component, and 
when instabilities occur, only the stellar mass necessary to 
restore stability is transferred from the disc to the bulge. 
No cold gas component is associated with the bulge in this 
model. In MORGANA , a significant fraction (half) of the 
baryonic mass (both gas and stars) of the disc is transferred 
to the bulge. This particular treatment avoids a series of 
consecutive instability events, that are instead frequent in 
the WDL08 model, in particular at high redshift. As we will 
show below, however, this modelling translates into a much 
more prominent role of the disc instability channel in bulge 
formation. 

We stress that both models adopted for disc instabilities 
are oversimplified, and provide a very crude description of 
the complex phenomenology of bar formation and evolution. 
In particular, the WDL08 model neglects the possibility that 
bar formation produces an inflow of gas towards the centre 
that could fuel starburst/AGN activity, and that can eventu- 
ally lead to bar disruption. The assumption of Vdisc = Knax, 
as well as the use of only the stellar mass disk in Eq. [T] are 
questionable. On the other hand, MORGANA makes more 
realistic assumptions about the disc circular velocity and 
includes the gaseous mass present in the disc in Eq. [1] In 
both models, the outcome of an instability event is mod- 
elled in a rather arbitrary way. We note that present sim- 
ulations do not provide clear indications about the fraction 
of disk mass that gets re-distributed, and how this depends 
on the halo/galaxy properties. In addition, the very crite- 
rion adopted to t ag a disk as unstable has been questioned 
in recent studies l| Athanassoulal [20081 ) . As we will show in 
the following, disk instability has important consequences 
on model predictions, and more work is certainly needed in 



order to improve this aspect of our modelling. In order to 
quantify the importance of this process, in the following we 
will also show model predictions obtained when the disc in- 
stability channel is switched off. We will refer to these runs 
as the pure mergers model runs. 

Another difference between the two models used in this 
study is given by how gas is treated during mergers. In the 
WDL08 model, the merger triggers a burst that converts a 
fraction of the combined gas into stars. These stars are added 
to the disc component of the remnant galaxy. In the MOR- 
GANA model, all gas and stars of the secondary are trans- 
ferred to the bulge of the remnant galaxy, and the cold gas 
associated with the bulge is efficiently converted into stars. 

Furthermore, we note that the WDL08 model accounts 
for satellite-satellite mergers, while MORGANA only consid- 
ers mergers between satellites and central galaxies. Finally, 
the small differences in the cosmological parameters adopted 
in the two models have little impact on model predictions. 

2.4 The Hopkins et al. prescriptions 

iHopkins et alj (2009a, HOP09 hereafter) analyzed a suite 
of hydrodynamic merger simulations and derived a 'gas- 
fraction dependent merger model'. We refer to the original 
paper for a detailed derivation of the model, while a sum- 
mary of_th^_Jcey_prescriptions can be found in Appendix 
A of IHopkins et all (2009b). In this model, the fraction of 
cold gas that participates in the starburst associated with a 
merger can be written as: 

^Tlburst 

burst -L 



m co M 



(1 + r crit /Rd) ■ exp(-r cri t/7?d) 



/ gas )-/disc-F(0,^-G(/x) 



where 

— = o-(l- 

ltd 

and /g as = m co id/(jn co id + m*,di 8 c) is the gas disc fraction, 
/disc = (m C oid + m„,disc)/mbar is the disc mass fraction, m ba r 
is the baryonic mass of the galaxy. Rd is the disc scale- 
length, and 6 is the inclination of the orbit relative to the 
disc. Assuming that, before coalescence, the distance of peri- 
centric passage is b = 2-Rd (typical of cosmological mergers), 
one obtains: 

fl ^ 

aF(0,u) = 

v ^ 1- 0.42 y/TTJi cosd 

where the parameter a subsumes details of the stellar profile 
shape and bar driven distortion dynamics during a merger. 
Finally, G(/x) contains the dependence on the merger mass 
ratio, and has the form: 

G(m) = -rj— 
1 + A* 

where n = 7712 /mi. In the literature, and even in the two 
papers by Hopkins et al. mentioned above, there are incon- 
sistencies in the definition of 'mass-ratio'. In the implemen- 
tation of the 'Hopkins' prescriptions used in this study, we 
define the mass of interest as the baryonic plus tightly bound 
central dark matter. Specifically, we have included the dark 
matter contribution in the following way: for each galaxy, 
we store the virial mass of the parent halo at the time of ac- 
cretion (i.e. at the last time the galaxy is central), and add 
10 per cent of this mass to the baryonic mass of the merg- 
ing galaxies. We have verified, however, that the inclusion 
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of the dark matter contribution does not affect significantly 
the results discussed below. The new stars formed during the 
merger driven starburst are added to the bulge component 
of the remnant galaxy. 

The whole stellar mass of the secondary is added to 
the spheroidal component of the remnant galaxy. It is fur- 
ther assumed that a fraction of the primary's stellar disc is 
transferred to the bulge component of the remnant galaxy, 
and is violently relaxed. Specifically, the mass of the disc 
that is 'destroyed' is: 

^^disc, destroyed /-^^'*,disc (2) 

We note that in the two standard models used in this study, 
the stellar disc of the primary is always unaffected during 
minor mergers and completely destroyed during major merg- 
ers. Therefore, we will monitor separately this model com- 
ponent when analysing results from the runs adopting these 
pres criptions. 

lliopkins et al.l (2009b) have investigated the implica- 
tions of the proposed model in a cosmological framework, 
using both empirical halo-o ccupation models and t he semi- 
analytic model presented in lSomerville et al.l (|2008l ). In par- 
ticular, they claim that their model leads to a significant 
suppression of bulge formation in low-mass galaxies, and 
that simulations and models that ignore the gas dependence 
of merger induced starbursts, have difficulties in reproducing 
the strong observed morphology-mass relation. 

We have implemented the prescriptions illustrated 
above in both models used in this study, and will re-address 
these issues below. We note that in the runs adopting the 
HOP09 prescriptions, we have always switched off the disc 
instability channel. In addition, in the WDL08 model, the 
runs using these prescriptions assume that the stars formed 
during merger-driven starbursts are added to the bu lge com- 
ponent of the remnant galaxy (as in lHopkins et al] 2009b). 
All other model details and parameters have been left un- 
changed. 



3 DEPENDENCY ON STELLAR AND HALO 
MASS 

In this section, we discuss the basic trends predicted by 
the two models used in this study, and analyze how they 
are modified by switching off the disc instability channel or 
adopting the HOP09 prescriptions discussed in the previous 
section. 

Fig. Q] shows the median (stellar) bulge-to-total ratio as 
a function of the galaxy stellar mass, at four different red- 
shifts. The left and middle columns show predictions from 
the WDL08 and MORGANA models, respectively, with dif- 
ferent colours used for different physical assumptions. The 
right column compares predictions from the two models. 

The WDL08 model (solid black lines) predicts a strong 
increase of the bulge-to-total ratio as a function of the galaxy 
stellar mass. When the disc instability channel for bulge for- 
mation is switched off (solid blue lines), the median bulge- 
to-total ratio decreases for all galaxies but the most massive 
ones which are still dominated by the bulge. Compared to 



this pure mergers run, the HOP09 on^f] predicts a larger me- 
dian bulge-to-total ratio for intermediate mass galaxies, but 
a smaller one for the most massive galaxies. In the WDL08 
model, the median bulge-to-total at fixed stellar mass de- 
creases slightly with increasing redshift. The scatter is large, 
as indicated by the shaded regions. This scatter is some- 
what reduced when considering central galaxies only, and it 
reflects the variation in galaxy (and halo) merger trees at 
fixed galaxy stellar mass. 

Also the standard MORGANA model (dashed black 
lines) predicts an increase of the median bulge-to-total ratio 
as a function of stellar mass, but this is somewhat shal- 
lower than that predicted by WDL08 . In the pure mergers 
run (blue dashed lines), the median bulge-to-total ratio at 
fixed stellar mass decreases with respect to the standard 
run, but it is higher than in the corresponding run of the 
WDL08 model, particularly at high redshift. The HOP09 pre- 
scriptions provide predictions that are very close to those of 
the pure mergers model for intermediate mass galaxies, but 
again lower bulge-to-total ratios for the most massive galax- 
ies. In the standard MORGANA run, the median bulge-to- 
total ratio increases slightly with increasing redshift, but 
with significant scatter at fixed stellar mass, as in WDL08 . 

The reason for the different behaviour obtained for in- 
termediate mass galaxies when adopting the WDL08 version 
of the HOP09 prescriptions can be ascribed to the different 
treatment of merger driven starbursts in the WDL08 model: 
during minor mergers, new stars are added to the disc com- 
ponent of the remnant galaxy in the standard run, while to 
the bulge component in the HOP09 run (see Section T2.3p . In 
contrast, in both the HOP09 and the pure mergers run of 
the MORGANA model, newly formed stars are added to the 
bulge component of the central galaxy. 

Interestingly, the standard MORGANA run predicts a 
quite large bulge-to-total ratio for galaxies of all masses, 
even at z ~ 2, where a very large fraction of the galaxies have 
B/T > 0.4. Even in the pure merger run, most galaxies have 
B/T > 0.2 at this redshift, and the median bulge-to-total 
ratio is significantly higher than for WDL08 (compare blue 
dashed and solid lines in the bottom right panel of Fig. [T]). 
Clearly, the efficient production of bulges in MORGANA at 
high redshift is not simply due to different assumptions made 
when the instability criterion is met, and is likely related 
to th e shorter merger time scales adopted ()De Lucia et al.l 
l201Cf ). We will come back to this issue later. 

Fig. [2] shows the median bulge-to-total ratio as a func- 
tion of the halo mass for the same redshifts as Fig.[T] The left 
and middle columns show results for each model only con- 
sidering central galaxies, while the right column compares 
all results for satellite galaxies only. 

Both models predict that a large fraction of central 
galaxies in haloes more massive than log[Mh a io] ^ 12.5 are 
dominated by the bulge component (the fraction of bulge- 
dominated central galaxies is significantly larger for MOR- 
GANA than WDL08 ). In both models, the HOP09 prescrip- 
tions provide a somewhat weaker increase of the bulge-to- 
total ratio as a function of halo mass, and a reduction of the 
median bulge-to-total ratio for the most massive galaxies. 



2 We recall that the disc instability channel has been switched 
off in the run adopting the HOP09 prescriptions. 
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Figure 1. Median bulge-to-total ratio as a function of the galaxy stellar mass. Different rows correspond to different redshifts, while 
different columns correspond to different models: the left column shows results from the WDL08 model; the middle column shows 
results from MORGANA , and the right column compares results from the two models. In all panels, black lines correspond to the 
standard models, blue lines correspond to the pure merger variant of these models, and red lines show results obtained using the HOP09 
prescriptions. The shaded regions in the left and middle columns show the 15th and 85th percentiles of the distributions obtained for 
the standard WDL08 and MORGANA runs respectively. 



The scatter at fixed halo mass is large, reflecting significant 
variations in the accretion histories of haloes (and of their 
central galaxies) at fixed mass. 



The median bulge-to-total ratio of satellite galaxies is 
relatively low (lower in the WDL08 model than in MOR- 
GANA ), and approximately constant as a function of halo 
mass. We note that, in these models, once a galaxy is ac- 
creted onto a larger system (i.e. becomes a satellite galaxy) 
its bulge-to-total ratio is unaffected, unless it suffers a 



merger with another satellite galaxj|f|. In the real Universe, 
tidal stripping and interactions with other satellite galaxies 
(e.g. harassment) are likely to increase the bulge-to-total ra- 
tio of satellite galaxies, increasing the median values plotted 
in Fig. [2] 

Interestingly, almost all central galaxies in haloes with 
mass slightly larger than ~ 10 12 M© are practically 'pure 
bulges' in the MORGANA model, and the vast majority 



3 Mergers between satellites are included in WDL08 but not in 
MORGANA . 
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Figure 2. Median bulge-to-total ratio as a function of the parent halo mass. Different rows correspond to different redshifts, while 
different columns and colours correspond to different models, as in Fig. [l] The left and middle columns have been obtained for central 
galaxies only, while the right column compares results for satellite galaxies. The shaded regions in the left and middle columns show the 
f 5th and 85th percentiles of the distribution for the standard and HOP09 runs. 



of central galaxies in 'Milky- Way type haloes' (with mass 
~ 1O 12 M ) have B/T > 0.6, at all redshifts considered. 
These results suggest that the standard MORGANA run has 
difficulties in forming a Milky- Way like galaxy in haloes of 
mass similar to that of our Galaxy Q. This does not ap- 
pear to be a problem in the standard WDL08 run (see also 
iDe Lucia fc Helmill200&1 ; iLi et al.|[20lb1 ). We note, however, 
that haloes of this mass are only marginally resolved in the 



4 iMaccio et"ai] fcoich have studied predictions from MORGANA 
for the luminosity function of 'Milky- Way' satellites. However, in 
the version of the model they use, the disc instability channel is 
switched off. 



simulations used in this study (with ~ 1000 particles in the 
WDL08 simulations). 



4 HOW AND WHEN DO BULGES FORM? 

The models we have in hands allow us to ask a number ques- 
tions about the formation of bulges: when did bulges form? 
Was most of their mass assembled during major or minor 
mergers? What is the relative importance of disc instabil- 
ity? How does this vary as a function of redshift? And in 
which environments did bulges form? 

In order to answer these questions, we have rerun our 
models and, each time the mass of the bulge is updated, 
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we have stored the information about the redshift, the halo 
mass and the fraction of mass contributed to the final bulge, 
distinguishing between major mergers, minor mergers, and 
disc instability. For the runs that use the HOP09 prescrip- 
tions, where we do not include the disc instability channel, 
we have stored separately the information about the bulge 
mass contributed through destruction of the primary's stel- 
lar disc (e.g. Eq. [2] in Section |2.4[) . We note that our def- 
inition of bulge formation refers to the event adding stars 
to the bulge of the selected galaxy or its main progenitor, 
i.e. the clock is reset for stars in a secondary galaxy once it 
merges with a more massive one. 

The 'channels' defined above correspond generally to 
a combination of different physical processes: e.g. in MOR- 
GANA , disc instability triggers both a re-arrangement of 
the stellar material originally distributed in the disc, and 
the inflow of disc gas towards the centre. This leads, in 
turn, to the formation of in situ new stars - a process that 
is not included in the WDL08 treatment of disc instability. 
Analogously, mergers are generally associated with both a 
starburst, and a re-arrangement of stars belonging to the 
merging galaxies. 

Fig. [3] shows the fraction of bulge mass contributed 
through different channels (different columns) for galaxies 
selected at z — as those having B/T > 0.4 (qualita- 
tively, the results do not change significantly when consid- 
ering all galaxies). We have split model galaxies into four 
stellar mass bins, and show the corresponding results in dif- 
ferent rows, with the most massive galaxies shown at the 
bottom. As in previous figures, solid lines show results from 
the WDL08 runs, while dashed lines show the corresponding 
results from MORGANA , with different colours referring to 
different runs. Where a black line is not visible, it overlaps 
perfectly with the corresponding blue line, meaning that 
switching off disc instability does not affect the contribu- 
tion of that particular channel to the final bulge mass (bot- 
tom panels). In the right column, the black lines show the 
fraction of bulge mass formed through the disc instability 
channel in the standard models, while red lines correspond 
to the fraction of bulge mass formed through destruction of 
the primary's disc in the runs that adopt the HOP09 pre- 
scriptions. We stress that these two sets of lines have been 
plotted in the same panels for convenience, and should not 
be compared against each other. 

Fig. [3] shows that the contribution from major merg- 
ers decreases with increasing stellar mass, while the con- 
tribution from minor mergers increases. The contribution 
from disc instability is largest for intermediate-mass galax- 
ies, and negligible for the most massive galaxies considered 
(bottom right panel). Both models and all runs considered 
share these trends, with a few notable differences: (i) bulges 
seem to form earlier in MORGANA than in WDL08 ; (ii) 
the contribution from disc instability is larger in the MOR- 
GANA model. For the most massive galaxies, disc instability 
contributes less than ~ 1 per cent of the final bulge mass, 
and all instabilities occurred at high redshift, in the small 
galaxies that later merged to form these massive systems. 
In MORGANA , the contribution from the disc instability 
channel is ~ 3 per cent for the most massive galaxies, but 
it tends to increase (weakly ) since z ~ 4 . Thes e results can 
be compared to those that iParry et al.l (120091 ) find for the 
semi-analytic model discussed in lBower et alj (|2006l ). Their 



Figure 8 shows that instabilities contribute to the bulges 
of present day galaxies significantly more than minor and 
major mergers, but for the most massive galaxies where 
the major merger contribution is dominant. Interestingly, 
they show that the contribution from disk instabilities in the 
'Durham' model is much larger than that in the 'Munich' 
model (which corresponds to the WDL08 model used in this 
study). The large contribution f rom disk instabilities in the 
Durham model is noted also in iBenson fc Devereuxl (|2010l . 
see e.g. their Figure 5), and is due to the assumption that in- 
stabilities result in the complete collapse of the disc (both of 
its stellar and gaseous component - see original paper for de- 
tails) . Clearly, the different outcome assumed for instability 
events has important consequences on the relative impor- 
tance of different channels to bulge for mation - a conclusion 
that s eems to contradict what found bv lBenson fc Devereuxl 
|20ld . appendix A2). We will come back to this issue later. 

In MORGANA , the HOP09 prescriptions result in a con- 
tribution from major and minor mergers that is approx- 
imately equal to that found in the standard run. In the 
WDL08 model instead, the HOP09 prescriptions result in a 
systematically lower contribution from major mergers, and 
higher contribution from minor mergers, at all redshifts. 
This is largely due to the fact that, when adopting the 
HOP09 prescriptions, the WDL08 model assumes that the 
stars formed during all merger driven starbursts are added 
to the bulge component of the remnant galaxy (rather than 
to the disc in the case of minor mergers, as in the stan- 
dard run). The reduced contribution from major mergers 
relates partly to the increased efficiency of bulge formation 
via minor mergers. One notable consequence of the HOP09 
prescriptions for the WDL08 model is for the most massive 
galaxies considered: for these, the major mergers channel 
becomes important only at z < 1. Finally, in both models, 
the disc destruction channel contributes significantly to the 
final bulge mass for the intermediate mass bins considered, 
but less than 10 per cent for the most massive galaxies in 
the sample. 

As mentioned above, the merger model adopted in 
MORGANA provides galaxy merger times th at are systemat- 
ically lower than those used in the WDL08 (|De Lucia et al.l 
l201Cf ). Since a significant fraction of the final bulge mass is 
associated with galaxy mergers (both minor and major), sys- 
tematic differences between the galaxy merger time-scales 
are expected to lead to a systematic difference in the char- 
acteristic formation times of galactic bulges. In order to un- 
derstand how these differences affect the results discussed 
above, we have re-run the MORGANA models using the 
same dynamical friction timescale prescriptions adopted in 
WDL08 . Results are shown in Fig. [4] We note that resid- 
ual merger times are assigned at the time of halo mergers 
in MORGANA , and orbital parameters are re-assigned after 
each major merger. In WDL08 , residual merger times for 
satellite galaxies are instead assigned when the parent dark 
matter substructures are stripped below the resolution limit 
of the simulation. So, although we are using now the same 
formulation of dynamical friction in the two models, overall 
merger times will still be different. 

The figure shows that, when using longer merger times 
for the MORGANA runs, bulges form later, particularly for 
the most massive galaxies considered. Interestingly, mak- 
ing merger time-scales longer also increases the contribu- 
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Figure 3. Fraction of bulge mass formed as a function of redshift through different channels (different columns), for galaxies with 
B/T > 0.4 selected at z = 0. Different rows correspond to different stellar mass bins, while different colours and linestyle correspond 
to results from different models as in Fig. ffj Where a black line is not visible, it overlaps perfectly with the corresponding blue line. In 
the right column, the black lines indicate the fraction of bulge mass formed through disc instability in the standard run, while red lines 
correspond to the fraction of bulge mass formed through destruction of the primary's stellar disc in the HOP09 prescriptions (see Eq. [2] 
in Section foil. 



tion from disc instability in the standard MORGANA run. 
This happens because galaxy discs now have longer times to 
develop instabilities. Increasing the merger times in MOR- 
GANA does not account for all differences between the two 
models used in this study. In particular, bulge formation still 



occurs earlier in MORGANA than in WDL08 for intermediate 
to low mass galaxies, and disc instability plays a much more 
important role in MORGANA than in the WDL08 model, 
particularly at high redshift. 
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Figure 4. As in Fig. [3] but using longer merger times in the MORGANA runs (see text for details). 



5 HOW AND WHERE DO BULGES FORM? 

In the previous section, we have analyzed when bulge for- 
mation occurs, and what is the relative importance of differ- 
ent channels at different times. Another question that can 
be addressed with our models is: what is the typical envi- 
ronment of bulge formation? Does it occur in groups or in 
the 'field'? And how does the characteristic environment of 
bulge formation vary as a function of cosmic time? 

In Figs \5\ and [6] we show the fraction of bulge mass 



contributed through different channels, as a function of red- 
shift and parent halo mass, for the WDL08 and MORGANA 
models respectively. Different rows correspond to different 
present day stellar mass bins, while different columns re- 
fer to different channels, as indicated by the legend. Data 
shown in Figs [5] and [6] have been computed for the stan- 
dard model considering all galaxies with B/T > 0.4, and 
have been normalized to the total bulge mass in each mass 
bin. Therefore, darker regions in each panel of Figs [5] and 
[6] indicate the ranges of redshift and halo mass where that 
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Figure 5. Fraction of the total bulge mass of galaxies of different stellar mass (increasing from top to bottom row) contributed from 
different channels (different columns), as a function of redshift and parent halo mass. For the standard WDL08 run. 



particular channel is more important. Qualitatively, the re- 
sults shown do not change when including all galaxies (i.e. 
without any cut for the bulge-to-total ratio). 

The figures show that the typical halo mass where 
different processes contribute to bulge formation increases 



with increasing stellar mass. Interestingly, for galaxies with 
log[M s tar] ~ 9 — 10, much of the bulge formation occurs 
in haloes of log[Mhalo] < 11-5 that are only marginally re- 
solved in the simulations used in this study. One clear differ- 
ence between the two models is that bulge formation occurs 
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Figure 6. Same as in Fig. [5] but for the standard MORGANA run. 



on a wider range of halo masses (more extended towards 
larger masses) in WDL08 . This difference is primarily seen 
for satellite galaxies. Indeed, when excluding these galax- 
ies, Figs [5] and Fig. [6] become more similar. For the major 
and minor mergers channels, this difference is due to the 
fact that the WDL08 model takes into account mergers be- 



tween satellite galaxies, that are not included in MORGANA 
. The contribution from satellite galaxies to the disc insta- 
bility channel comes from the fact that in WDL08 , the disc 
radius of a satellite decreases in proportion to the radius of 
its dark matter halo. This might not be generally true and 
could artificially increase the bulge-to-total ratio of satellite 
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galaxies. We have verified that fixing the disk radius at the 
time of infal{3 decreases significantly the contribution from 
disk instability due to satellite galaxies. We note, however, 
that our model does not include physical processes such as 
tidal stripping and harassment that would again increase 
the bulge-to-total ratio of satellites. 

Another obvious difference between Figs [5] and [5] is the 
more pronounced (and extended) contribution from disc in- 
stability to the formation of the most massive bulges in the 
standard MORGANA run with respect to the WDL08 model 
(see also bottom right panel of Fig. [3]). Finally, these figures 
show that there is a somewhat 'tighter' correlation between 
the mass of the halo and the redshift in MORGANA . This is 
likely due to the fact that mass accretion histories obtained 
using PINOCCHIO are 'smoother' than those obtained from 
numerical simulations. 



6 ELLIPTICALS AND DISC REGROWTH 

In this section, we will focus on galaxies that are dominated 
by a bulge, and that we will call 'ellipticals' (for a more 
extended discussi on on the formation history of elliptical 
galaxies, see also |Pe Lucia et"afl l2006l fl. More specifically, 
we include in the elliptical class all galaxies with at most ten 
per cent of the stellar mass in a disc (B/T > 0.90). We will 
address, in particular, three specific questions: (i) what is the 
typical stellar mass and environment of elliptical galaxies? 
(ii) what is the frequency and relevance of disc regrowth for 
these galaxies? (iii) when do galaxies become ellipticals and 
through which physical process(es)? 

Fig. [7] shows the fraction of galaxies classified as ellipti- 
cals, as a function of the galaxy stellar mass (left panels) and 
of the parent halo mass (right panels) . Top and bottom pan- 
els refer to the WDL08 and MORGANA models, respectively. 
As the stellar mass increases, a larger fraction of galax- 
ies are classified as ellipticals, and the distributions com- 
puted from MORGANA extend to larger masses than those 
obtained for WDL08 . We note that an excess of massive 
galaxies in the MORGANA model is well documented and is 
primarily due to an inefficient suppression o f star formation 
in central galaxies by AGN feedback (e.g. iFontanot et al.l 
12009] ; li< imm et alJ [2009). Disc instability does not affect sig- 
nificantly the number and distribution of elliptical galax- 
ies in WDL08 , while the fraction of ellipticals is sensibly 
reduced when switching off the disc instability channel in 
MORGANA , also at relatively large masses. The distribu- 
tion as a function of halo mass is approximately flat for the 
WDL08 model, and more skewed towards low-mass haloes 
in MORGANA (i.e. a larger fraction of galaxies residing in 
relatively low-mass haloes are classified as ellipticals in this 
model), where a large fraction of ellipticals in relatively low- 
mass haloes form through disc instability. The decreasing 
fraction to higher halo mass in this model is due to the 
increasing contribution of disc-dominated satellite galaxies. 



5 IWeinmann et al. have verified that results from the 
model would be virtually unchanged but for the morphology of 
satellite galaxies. 

6 Note, however, that a different definition of 'ellipticals' was 
adopted in that study. 
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Figure 8. Top panel: Mean regrowth rate (see text for details). 
Bottom panel: fraction of central galaxies that regrow a disc. 
Black lines refer to the standard WDL08 run, while red lines 
correspond to the standard MORGANA run. Different linestyles 
correspond to galaxies of different stellar mass, as indicated in 
the legend. 



Interestingly, the HOP09 prescriptions tend to reduce signif- 
icantly the number of ellipticals in both models. Overall, the 
number densities of galaxies that are classified as ellipticals 
in MORGANA are much larger (more than a factor two) than 
those obtained in WDL08 . We note that model results are 
in qualitative agreement with observational data indicating 
that the total fraction of ellip ticals is not expecte d to be 
higher in more massive haloes l|Wilman et alj|2009h . 

As discussed in Sections I2.ll and I2.2I bulge dominated 
galaxies can regrow a new disc, provided there is enough 
cooling from the surrounding hot halo gas. In principle, mi- 
nor mergers with gas-rich satellites also form new disc stars 
in the WDL08 model. However, most of the galaxies accreted 
onto centrals of relatively massive haloe s (where most of the 
centr al ellipticals are) will be gas-poor (|De Lucia fc BlaizotJ 
l2007f ). and mergers between satellites are rare. Therefore, 
the only mechanism through which galaxies can grow a new 
disc is by accretion of fresh gas material from the surround- 
ing hot halo. 

In order to quantify the importance of disc regrowth, we 
have selected all central elliptical galaxies from our models 
and analyzed their merger trees, storing the increase of their 
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Figure 7. Distribution of galaxies with B/T > 0.90 at z = 0, as a function of the galaxy stellar mass (left panels) and of the parent 
halo mass (right panels). The top and bottom rows show results for the WDL08 and MORGANA models respectively. Different colours 
correspond to different runs: black lines correspond to the standard models, blue lines correspond to the pure mergers variant of these 
models, and red lines show results obtained using the HOP09 prescriptions. The number densities of galaxies are given in the legend of 
the right panels, and are expressed in units of Mpc - 3 . The histograms have been normalized dividing by the number of galaxies in each 
stellar or halo mass bin. 



stellar disc mass as a function of time, until they become 
satellites (or until z — 0, for central galaxies). The top panel 
of Fig.[8]shows the mean regrowth rate for galaxies in differ- 
ent stellar mass bins (the stellar mass corresponds to that 
of the galaxy at the redshift under consideration). The rates 
shown in Fig. [8]have been normalised to the total number of 
(central) galaxies experiencing regrowth. The figure shows 
that, on average, the rate of disc regrowth decreases with 
decreasing redshift, for galaxies of all masses. For the most 
massive galaxies considered (dot-dashed lines in Fig. [5]), the 
rate of regrowth averaged over 1 Gyr time-scale is always sig- 
nificantly smaller than the galaxy mass at the corresponding 
time. This means that a galaxy that crosses the threshold 
B/T — 0.90 when it is already rather massive, will likely 
stay above this threshold at any later time. For lower mass 
galaxies, the rates are more significant, particularly at high 
redshift where galaxies are more gas-rich. On average, how- 
ever, these galaxies will not grow a large disc: a galaxy with 
stellar mass ~ 10 10 Mq and B/T = 1 will regrow a disc 
containing about 10 per cent of the mass in ~ 3 — 7 Gyr 
at z ~ 0, or ~ 1 — 3 Gyr at z ~ 1. Disc regrowth is also 
inhibited by subsequent mergers which increase the bulge- 
to-total ratio. These results suggest that disc regrowth is not 
significant for the most massive galaxies, but that it repre- 



sents a non negligible component in the evolution of low and 
intermediate mass galaxies, particularly at high redshift. 

The bottom panel of Fig. [S] shows the fraction of cen- 
tral galaxies experiencing regrowth in the two models used 
in this study. At the highest redshifts considered, the bal- 
ance between cooling and feedback is such that more galaxies 
grow a disc in WDL08 than in MORGANA . As the redshift 
decreases, AGN feedback becomes more and more impor- 
tant, particularly for the most massive galaxies (that are 
sitting in the most massive haloes). This determines a sig- 
nificant decrease of the fraction of central galaxies experi- 
encing regrowth in WDL08 : only about 20 per cent of the 
most massive galaxies are growing a new disc, and, as ex- 
plained above, this occurs at relatively low rates. On the 
contrary, the fraction of central galaxies experiencing disc 
regrowth increases with decreasing redshift in MORGANA . 
Albeit the regrowth rates are low also in this model, almost 
all galaxies (~ 90 per cent of the most massive ones) are 
growing a new disc. This different behaviour can be ascribed 
to the different treatment of radio-mode AGN feedback, that 
is much more effici ent in suppressing cooling in WDL08 tha n 
in MORGANA (e.g. lKimm et al.ll2009l ; lFontanot et al.ll2010l ). 

Another interesting question to address about elliptical 
galaxies is when and through which mechanism(s) they ac- 
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Figure 9. Distribution of the times when galaxies cross the threshold B/T = 0.9 for the first time. Galaxies are split by their final 
stellar mass (different rows) and each column shows the contribution from different channels. Galaxy counts are weighted by the fraction 
of bulge mass contributed by each process, and histograms are normalised by the total number of galaxies in each mass bin. Black and 
red lines are used for the WDL08 and MORGANA model, respectively. Red dashed lines show results obtained from MORGANA using 
longer merger times. 



quired their morphology. We address this question in Fig. [9] 
where we show the redshift at which galaxies cross the 
threshold B/T = 0.9 for the first time. Galaxies are split by 
their final mass and different columns correspond to differ- 
ent bulge formation channels. When several processes con- 
tribute to make the galaxy cross the threshold B/T = 0.9, 
we have weighted the counts by the fraction of bulge mass 
contributed by each channel. The figure shows that most 
galaxies acquire an 'elliptical' morphology because of a ma- 
jor merger event (left column), while minor mergers seem 
to play a negligible role in turning galaxies into ellipticals, 
particularly in WDL08 . Disc instability is responsible for 
turning a relatively large fraction of galaxies into ellipticals 
in MORGANA , particularly at intermediate masses, while 
it plays a much less prominent role in WDL08 . The distri- 
butions obtained from MORGANA are peaked at redshifts 
significantly higher than in WDL08 , as a consequence of the 



significantly shorter merger times and of the more promi- 
nent role of disc instability. Dotted lines in Fig. [9] show 
results from MORGANA obtained adopting longer merger 
times (see Section [4}, and confirms that galaxies become 
ellipticals later when longer merger times are adopted. As 
discussed in Section 3] the figure also shows that this run 
is characterized by a larger contribution from disc instabil- 
ity to bulge formation. We stress that the times plotted in 
Fig. [5] should not be confused with the 'formation times' of 
elliptical galaxies, as these times are not related to the star 
formation history of these galaxies. 

Not all galaxies considered in Fig. [5] are still ellipticals 
at z = 0. As discussed above, however, the regrowth rates 
are lower for intermediate-massive galaxies so that a larger 
fraction of the most massive galaxies maintain their elliptical 
morphology down to z = 0. In particular, we find that in the 
WDL08 model about 30, 64, and 96 per cent of the galaxies 
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in each of the mass bins considered (in order of increasing 
mass) still have B/T > 0.9 at redshift zero. In MORGANA , 
the corresponding fractions are 85, 46, and 88 per cent. The 
fractions increase to about 53, 79, and 98 for the WDL08 
model and 89, 61, and 92 per cent for MORGANA , when 
considering galaxies with B/T > 0.6 at z = 0. 



7 DISCUSSION AND CONCLUSIONS 

In this paper, we have analyzed predictions for the 
formation of bulges from two independently developed 
galaxy formation models. In particular, we have con- 
sidered (i) the rece nt implementation of th e Munich 
model (WDL08 ) by |Pe Lucia fc Blaizotl l|2007f ). with its 

generalizatio n to the WMAP3 cosmology discussed in 

I Wang et all (l2008f). and (ii) the MORGA NA model presented 

to a 



Monaco. Fontanot fc Taffonil J2007l). and adapted 



WMAP3 cosmology as described in lLo Faro et ail (|2009T ). 

The two models include the same channels for bulge 
formation (mergers and disc instability), but assume dif- 
ferent prescriptions to model these physical processes. In 
order to study how results vary as a function of specific 
physical assumptions, we have also implem ented alternative 
merger prescriptions l|Hopk ins et al. 2009), based on results 
from recent hydrodynamic merger simulations. In this pa- 
per, we have focused on theoretical predictions, postponing 
to a forthcoming paper a detailed comparisons with observa- 
tional results (Wilman et al., in preparation). In a compan- 
ion paper (Fontanot et al., in preparation), we will study the 
physical properties and formation histories of galaxies with 
no significant bulge component, as predicted by the same 
models considered here. 

Both models used in this study, with all different physi- 
cal assumptions considered, predict a strong correlation be- 
tween the galaxy morphology and its stellar mass, with more 
massive galaxies having larger bulge-to-total ratios. For cen- 
tral galaxies, there is also a strong correlation between the 
morphology of a galaxy and its parent halo mass, with 
most of the central galaxies of haloes with mass larger than 
~ 10 13 Mq being dominated by a bulge. These trends are 
not surprising, given our assumption that bulges form dur- 
ing mergers, and the strong correlation between the galaxy 
mass and the mass of the parent halo for central galaxies: 
more massive galaxies will generally have a richer merger 
history than their less massive counterparts, and more mas- 
sive galaxies will sit at the centre of more massive haloes. 

Taking advantage of our models, we have studied in de- 
tail the contribution to bulge formation from different 'chan- 
nels' (major and minor mergers, and disc instability). Dif- 
ferences arise between the different models and implemen- 
tations, but the results at redshift zero can be summarized 
as follows: 

(i) major mergers dominate the contribution to bulges of 
galaxies less massive than ~ 10 10 M o ; 

(ii) for galaxies more massive than ~ 10 10 Mq , the contri- 
bution from minor and major mergers are comparable; 

(iii) disc instability represents the dominant contribution 
to the formation of bulges of intermediate mass galaxies 
(~ 10 10 - 10 11 M ). 

Qualitatively, our results are in agreement with what 



found by iParrv et all (|2009h . although in their model disk 
instability plays a more prominent role. It is worrying that 
such an important contribution to bulge formation comes 
from the process that we probably model in the poorest 
way (disc instabilities). The results discussed in this paper 
confirm that further work is needed in this area in order 
to improve our galaxy formation models. This is true not 
just for the criterion ado pted to tag a disc a s unstable (as 
discussed for example in I Athanass oula 2008|), but also for 
the treatment of these events. In fact, the two models used 
in this study assume the same mathematical criterion for 
disc instability but adopt quite different assumptions for the 
physical quantities considered, and model the outcome of in- 
stability events in different ways. In particular, the WDL08 
model only transfers to the bulge a fraction of the stellar disc 
that is enough to restore stability. In MORGANA , half of the 
baryonic mass of the disc (i.e. gas and stars) is transferred 
to the bulge component each time an instability episode oc- 
curs. At high redshift, where galaxies are more gas rich and 
dynamical times are shorter, inflow of gas towards the centre 
leads to a rapid and efficient growth of the central spheroidal 
component in this model. 

It is instructive to see how different the adopted models 
are. We can estimate the fraction of discs that are unstable 
by looking at the distributions of the left-hand side of Eq. [1] 
for the pure merger models. Results are shown in the bottom 
panel of Fig. 1101 while the other two panels show the distri- 
butions of the disc velocity and radius entering Eq. [T]in the 
same models. The figure shows that there are systematic 
differences between these quantities in the two models. In 
particular, the disc velocities and disc scale-lengths used in 
MORGANA are systematically larger and smaller than those 
adopted in WDL08 , respectively. Despite these differences, 
the fractions of unstable discs in the two models are com- 
parable: integrating the distributions shown in the bottom 
panel up to en m = 0.75 for WDL08 and en m = 0.7 for MOR- 
GANA (see Section [2}, one obtains that ~ 5 — 7 per cent 
of the discs are unstable in both models. This fraction is 
approximately independent of redshift. The different treat- 
ment of instabilities, however, leads to a much more promi- 
nent role of the disc instability channel in MORGANA . To 
further test our conclusion, we have verified that the con- 
tribution from disc instability (particularly at high redshift) 
is significantly reduced in MORGANA when the fraction of 
disc mass transferred to the central spheroidal component is 
reduced from 0.5 to 0.1. 

In the framework of our models, bulge dominated galax- 
ies can grow a new disc if they are fed by appreciable cool- 
ing flows. The rates of disc regrowth are relatively low at 
low redshift and for massive galaxies. They are, however, 
more significant for intermediate and low mass galaxies and 
at high redshifts, where ellipticals can change back their 
morphology and become disc dominated galaxies at some 
later time (i.e. the morphology of these galaxies is transient). 
When an efficient radio-mode feedback is assumed (like in 
the WDL08 model used here), there is a decline of the typical 
regrowth rate of the most massive galaxies at low redshift, 
and of the fraction of these galaxies experiencing regrowth. 
This is expected given that these galaxies live in the most 
massive haloes where the radio-mode is assumed to play a 
major role. 

Only a minor fraction of the ellipticals in our models ac- 
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Figure 10. Distributions of the disc velocity and radius entering 
Eq. [l] at two different redshifts. Black lines correspond to the 
standard WDL08 model, while red lines are for MORGANA . 
The bottom panel shows the distribution of the left-hand side of 
Eq. [l] and the vertical lines mark the threshold for unstable discs 
adopted in the two models. 



quire their morphology through minor mergers. In WDL08 , 
the vast majority of the galaxies become ellipticals through 
major mergers, and this occurs at relatively low redshift 
(z < 2, only a very small fraction of the galaxies become 
ellipticals at higher redshift). In MORGANA , the change of 
morphology occurs at higher redshift, and a large fraction 
of the galaxies cross the threshold we have adopted to de- 
fine galaxies as ellipticals (B/T = 0.9) as a consequence of 
disc instabilities. We have demonstrated that this is due to 
a combination of significantly shorter merger times and a 
different treatment of disc instability events. 

The implementation of a gas-fraction dependent merger 
model provides trends in bulge-to-total ratio that do not de- 
viate strongly from those of our pure mergers runs at inter- 
mediate masses. This appears to be in con tradiction with 
the previous claim of iHopkins et all {2009) who find that 
taking into account the gas dependence of merger induced 
starbursts reduces bulge formation in galaxies less massive 



than ~ 10 10 Mq . We note, however, that the results from the 
'simplified model' discussed in Hopkins et al. are not repre- 
sentative of the results from our semi-analytic models. In 
particular, their simplified model provides much more mass 
in bulges at galaxy masses lower than ~ 10 10 M@, where the 
assumption of a gas dependent merger model should make 
the largest difference. We note that our standard models 
do include a dependency on the gas available during mi- 
nor mergers, and they do so in different ways. The WDL08 
model assumes that the stars formed during starbursts asso- 
ciated with these mergers are added to the disc component 
of the remnant galaxy, while MORGANA adds these stars to 
the bulge component. The results discussed above show that 
this model difference does not significantly influence model 
results, because minor mergers represent a relatively minor 
contribution to bulge formation. 

As explained above, our models naturally predict a 
strong correlation between the galaxy stellar mass and its 
morphology. In addition, this correlation evolves in such a 
way that fewer galaxies are bulge dominated at higher red- 
shift. One of the models used in this study (MORGANA ) 
does predict an 'excess' of galaxies with large bulge-to-total 
ratios at intermediate-low masses. As emphasized above, 
however, this is not due to the fact that the model neglects 
the dependency on gas fraction. Rather, the behaviour of 
MORGANA can be ascribed to galaxy merger times that 
are significantly shorter than those adopted in WDL08 , and 
to a different treatment during instability events that leads 
to an efficient production of bulges at high redshift. In- 
terestingly, the implementation of a gas-fraction dependent 
merger model reduces significantly the number of bulge dom- 
inated galaxies (ellipticals), and delays their formation. This 
is a consequence of the assumption that some fraction of the 
stellar disc is always preserved, even during major mergers. 

It remains to be seen if any of the models discussed 
in this paper provides a good agreement with the grow- 
ing amount of information being accumulated on the cos- 
mic mophological mix, at different epochs. Such a compar- 
ison will certainly represent an important test-bed for the 
next generation of models, and provide stringent constraints 
on the mechanisms through which bulges are assembled. 
While we defer to a forthcoming paper a detailed compar- 
ison with observational data, we note that the predicted 
bulge-to-total distributions (as a function of both stellar and 
halo mass - see Section [3]), as well as the total fraction of 
bulge-dominated galaxies, are affected by the physical pre- 
scriptions adopted to model bulge formation. None of these 
observables has been used to 'tune' the models in the first 
place, so that they can be considered as genuine model pre- 
dictions and tested against observational measurements. 

We have shown that the contribution to bulge mass 
from different channels is also affected by the adopted phys- 
ical modelling. If 'classical bulges' are primarily formed 
through mergers and 'pseudo-bulges' can be associated with 
secular evolution, the results discussed in this paper provide 
predictions for the relative importance of these two popu- 
lations, at different cosmic epochs and in different environ- 
ments. Distinguishing pseudo-bulges from classical bulges is 
difficult, requires good photometric data and, ideally, also 
high- resolution spectroscopic information (Deb attista et all 
2005). S ome st atistics are, however, available at low redshift 
l Gadottill2009T ) and can be used to constrain our models. 
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